Is the concept of the non-Hermitian effective Hamiltonian 
relevant in the case of potential scattering? 
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We examine the notion and properties of the non-Hermitian effective Hamiltonian of an unstable 
system using as an example potential resonance scattering with a fixed angular momentum. We 
present a consistent self-adjoint formulation of the problem of scattering on a finite-range potential, 
which is based on separation of the configuration space into two segments, internal and external. 
The scattering amplitude is expressed in terms of the resolvent of a non-Hermitian operator 7i. The 
explicit form of this operator depends on both the radius of separation and the boundary conditions 
at this place, which can be chosen in many different ways. We discuss this freedom and show 
explicitly that the physical scattering amplitude is, nevertheless, unique although not all choices are 
equally adequate from the physical point of view. 

The energy- dependent operator Ti should not be confused with the non-Hermitian effective Hamil- 
tonian TtcH which is usually exploited to describe interference of overlapping resonances. We note 
that the simple Breit-Wigner approximation is as a rule valid for any individual resonance in the case 
of few-channel scattering on a flat billiard-like cavity, leaving no room for non-trivial Tied to appear. 
The physics is appreciably richer in the case of an open chain of L connected similar cavities whose 
spectrum has a band structure. For a fixed band of L overlapping resonances, the smooth energy 
dependence of 7i can be ignored so that the constant L x L submatrix TL c s approximately describes 
the time evolution of the chain in the energy domain of the band and the complex eigenvalues of 
"W e ff define the energies and widths of the resonances. We apply the developed formalism to the 
problem of a chain of L 8 barriers whose solution is also found independently in a closed form. We 
construct 7i c ff for the two commonly considered types of the boundary conditions (Neumann and 
Dirichlet) for the internal motion. Although the final results are in perfect coincidence, somewhat 
different physical patterns arise of the trend of the system with growing openness. Formation in the 
outer well of a short-lived doorway state shifted in the energy is explicitly demonstrated together 
with the appearance of L — 1 long-lived states trapped in the inner part of the chain. 

PACS numbers: 05.60.Gg, 03.65.Nk, 24.30.-v, 73.23.-b 



I. INTRODUCTION 

To the best of our knowledge, the concept of the non- 
Hermitian effective Hamiltonian appeared first in Fesh- 
bach's papers 0, 01 m connection with the general the- 
ory of resonance nuclear reactions 0, Hand, indepen- 
dently, in Livsic's study of open systems 0, . A typical 
atomic nucleus forms near an excitation energy E above 
the threshold of nucleon emission a rich set of long-lived 
compound states with a very dense energy spectrum. For 
relatively small excitation energies, these states manifest 
themselves as narrow isolated resonances in collisions of 
nucleons and target nuclei. The sharp energy dependence 
of the corresponding cross sections near a given resonance 
is described by the universal Breit-Wigner formula. Any 
smoother variations can be ignored in the domain of an 
isolated resonance. At higher E, the resonance states be- 
gin to overlap and strongly interfere. Nevertheless, one 
can still neglect the smooth energy dependence within 
a group of close interfering resonances. In this approxi- 
mation, the propagation of the unstable system created 
on the intermediate stage of the collision is characterized 
by the resolvent of an energy-independent non-Hermitian 
matrix [2J. The latter describes in the time picture the 
irreversible evolution of the excited intermediate state 



and can therefore be interpreted as the effective Hamil- 
tonian 7i c ff whose anti-Hcrmitian part, originating from 
on-shell self-energy contributions, is responsible for de- 
cays into open channels. The complex eigenvalues of the 
effective Hamiltonian determine the energies and widths 
of the resonances. The method of the effective Hamilto- 
nian proved to be a success for describing many impor- 
tant resonance phenomena , espec i ally in the context of 
chaotic scattering @, ©, [H IH 111 III El 

It must be stressed that, similar to the Breit-Wigner 
formula, the notion of the effective Hamiltonian is local 
in energy although the general approaches |3|, la work 
in much wider energy intervals where smooth variations 
become important and should already be taken into ac- 
count. For the very concept of the effective Hamilto- 
nian to be consistent, the scales of resonant and smooth 
variations must be appreciably different. Otherwise the 
resonance behavior of scattering amplitudes will be dis- 
torted or even completely destroyed. The smooth depen- 
dence influences the background phases as well as the 
parameters of the resonances situated in different energy 
domains. Such effects cannot be described by simply 
enlarging the dimension of the matrix of the effective 
Hamiltonian. Instead, a large energy-dependent matrix 
Ti.(E) emerges whose simple interpretation as a time-shift 
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operator is no longer valid. There is, generally, no one- 
to-one correspondence between the complex energies of 
resonance states on the one hand and the independent 
eigenvalues of this large matrix on the other. 

The current intensive studies, theoretical as well as 
experimental, of chaotic scattering of a particle by an 
open two-dimensional cavity (see, e.g., jl5| and refer- 
ences therein) renewed interest in the Hamiltonian ap- 
proach to the resonance scattering theory and, in par- 
ticular, to the concept and properties of the effective 
non-Hcrmitian Hamiltonian [If* Il7l Il8| . Although the 
scattering processes considered are purely potential, the 
theory is formulated in close analogy with the formal- 
ism developed in Q, S for nuclear collisions. The 
configuration space is divided into two parts, internal and 
external. The Feshbach projection technique is employed 
to express the scattering matrix in the specific form that 
explicitly reveals the sharp resonant energy dependence 
while all smoother variations show up only indirectly via 
changes of the matrix elements. Related is the problem 
of electromagnetic field quantization in open optical cav- 
ities [l9|. 

While the probability amplitudes of the physical pro- 
cesses are fixed unambiguously, the projection procedure 
exploited is not unique. There exists a rather wide free- 
dom in choosing the surface of separation as well as the 
boundary conditions (BCs) on this surface. Therefore, 
the amplitudes of interest are expressed in terms of quan- 
tities that depend on the details of the formalism. The 
independence of the final results of the calculations is 
not, as a rule, directly seen. Therefore, a certain cau- 
tion is necessary to avoid incorrect assertions. An un- 
expected dependence on the (auxiliary) boundary con- 
dition at the cavity-lead interface was found, for exam- 
ple, in 17]. Moreover, poor agreement with numerically 
calculated exact poles of the S matrix was revealed for 
poles distant from the real axis. These points, partly 
attributed by the authors to numerical limitations, re- 
quire further clarification. Additional physical arguments 
should sometimes be involved to reasonably restrict the 
freedom. 

It is important to recognize that the density of levels in 
billiard- like cavities is much lower than that of the many- 
body nuclear systems. Actually, strong overlap and in- 
terference of different resonance states are not, as a rule, 
possible in the single two-dimensional cavities ordinarily 
considered 01- The simple Breit-Wigner approximation 
is usually sufficient in this case for any individual res- 
onance (see, e.g., [53, EJ)- Noticeable interference can 
occur only in the rare event of accidental near degeneracy 
of resonances. Experimental observation of such an in- 
terference of few resonances in an open microwave cavity 
has recently been reported in |22fl. Stronger overlap and 
a non trivial effective Hamiltonian matrix can, however, 
appear when open chains of a number of similar poten- 
tial wells connected to each other are considered. The 
energy spectrum has a band structure in this case with 
much denser spectrum within a given band. A schematic 



model of such a kind was considered in [23j. A similar 
example was also investigated in 0] m the framework of 
graph theory j24[ . Some general aspects of scattering in 
periodic structures have also been considered in |25j |. 

In this paper we examine the notion of the effective 
Hamiltonian in potential scattering. The questions we 
are concerned with do not depend on the regular or irreg- 
ular character of the motion. Thus, we restrict ourselves 
to the simplest case of the single- channel s-wave scatter- 
ing. The extension to higher partial waves is straightfor- 
ward. A consistent self-adjoint formulation of the prob- 
lem of scattering on a finite-size potential in terms of the 
internal and external subproblems is presented in Sec. II. 
The (one-dimensional) S matrix is expressed in terms 
of a non-Hermitian energy-dependent operator TL whose 
form, together with the form of the associated R func- 
tion (Sec. Ill), depends on the BC and the radius a of 
separation in the configuration space. Different BCs de- 
fine different representations of the internal and external 
parts of the (unique) scattering wave function. The role 
of the fictitious direct reflection at the separation point 
a is also discussed. 

In Sec. IV the exactly solvable open Kronig-Penney 
model is considered as an example of potential scatter- 
ing with interfering resonances. Closed expressions are 
found for both scattering S(k) and staying-wave K{k) 
functions as well as Wigner's function R(k). Their an- 
alytical properties in the complex k and energy planes 
are analyzed in detail. Hereafter we compare these exact 
results with those obtained in the framework of the pro- 
jection formalism of the previous sections. We analyze 
two typical choices of the BCs for the internal motion: 
the cases of Neumann and Dirichlet BCs. In the second 
case, the internal problem corresponds to a closed coun- 
terpart of the system under consideration that allows one 
to follow the changes of the motion due to the interference 
of resonances when the openness grows. The latter is in 
line with numerous applications considered in the liter- 
ature. Non-Hermitian effective Hamiltonians for a fixed 
band of resonance states are built up in both cases. We 
summarize our main findings in the concluding Sec. V. 



II. SEPARATION OF THE HILBERT SPACE 

The radial motion in the s-wave scattering is described 
by the Schrodinger equation (we use the units ti 2 /2m = 1 
throughout the paper) 

(^^2+U(r^x(r) = k\(r) (1) 

with the boundary condition x(0) = at the origin. We 
use below the specific method which, basically, goes back 
to Bloch's paper [2(| (see also [13 )• However, we present 
a derivation which leads to interrelated boundary prob- 
lems for internal and external regions. Making use of 
Heaviside's step function one can decompose the wave 
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function into internal and external parts as follows 
x(r) = u(r) 9(a — r) + 4>{r) 9{r — a) , 



(2) 



where the functions u and <f> are supposed to be continu- 
ous at the point a together with their two first derivatives. 
The partition radius a can be chosen arbitrarily; we sug- 
gest only that U(r > a) = 0. The second derivative of 
the function (J2J) reads 



X"(r)=u"(r), 



r) + <j)"{r) 6(r - a) + 



[4>'{a) - u'{a)} S(r - a) - [u(a) - 0(a)] S'(r - a) . (3) 

When substituting this expression into Eq. Q), there ex- 
ists a freedom in attributing the local terms in to the 
internal or external regions. One possible choice is to de- 
fine the singular functions as S(r — a_) = (r|a_) = (&_|r) 
and S'(r - a+) = S(r - o+) = -<r|a+>' = -{a+\'r), 
where a± = a ± 0, thus considering the first as belonging 
to the internal subspace and the second to the external 
one, respectively. 

With such definitions, we can represent Q in the form 



H, 



with the following entries: 




= k 2 



H, 



(N) 



H. 



K + U+\a-){a-\'. 
= k+\a,)'(a,\, 



V = -\a-){a+\' , Vi = -\a+y(a-\ 



(4) 



(5a) 

(5b) 
(5c) 



where K stands for the kinetic energy operator. In par- 
ticular, in the coordinate representation the matrix ele- 
ments of the entries take the form 



(r\V\r') 



■£x + U(r)\ 6(r-r') 
-5(r- a-)5'(r' - a_) , (6a) 

~6'(r~a+)6(r' - a+) , (6b) 
(r'|yt| r ) = 6 (r - a-) 5'{r' - a+) .(6c) 



The presence of the singular terms assures Hermiticity of 
the operators iJ; n t and H cy ± in contrast to K alone [26j. 
The relations (liil-Hint^a) = {H- m tUi\u2) , etc., can be eas- 
ily checked by means of partial integrations. The singular 
operators on the boundary provide the boundary condi- 
tions, the first one Neumann for the internal region, the 
second one Dirichlet for the external region. Equation © 
implies that both BCs are interrelated. This also leads 
to the adjointness of the coupling operators V and V*. 
The Dirichlet boundary operator is essentially different 
from that of Neumann type and cannot be produced by 
the boundary operators as used in pr| . 

The range of operators in Eq. JfJJl within the Hilbert 
space is defined by the functions on which the singular 



terms vanish. This requirement fixes the BCs at the point 
of separation a. The full Hilbert space is a direct sum 
of the space spanned by the eigenvectors of the internal 
problem < r < a with Neumann BC, 

(-^ + C/ m) <(r)=4 N) <M, <'(«) = (7) 

and that spanned by the eigenvectors of the external 
problem r > a with Dirichlet BC, 



dr 2 



0°(r) = P0°(r), 02(a) =0; fc>0. (8) 



The type of BC is explicitly indicated by the correspond- 
ing superscripts in Eqs. Q-lJBJ. As usual, we assume 
normalization conditions 



u° n {x)u° m {x) dx = S ml 



<t>l(x)4> k ,(x) dx = S(k-k') 



(9a) 
(9b) 



in the discrete and continuous spectra, respectively. The 
external function 



sm[k(r — a)] oc e 



-ikr 



-2ika Akr 



(10) 



describes a wave fully reflected at the point a, the corre- 
sponding S function being equal to Sq (k) = e~ 2lka . 

One can also proceed in the opposite way and ascribe 
the terms with 5' and with S to the internal and external 
regions, respectively. In this case, the entries in Eq. Q 
read 



V 



H^ t =K-\a + )(a+\', 
a-Y(a+\, V^ = \a+){a. 



(11a) 
(lib) 
(11c) 



That corresponds to the nonperturbed internal problem 
with Dirichlet BC u„(a) = 0, when the external one has 
Neumann BC tp® (a) = 0. The latter results in an addi- 
tional shift by 7r/2 of the reflection phase at the separa- 
tion point a: S^(k) = S ( D) (k) = e -2i(fca +3 r/2) . 

The interplay of the internal and external motions due 
to the off-diagonal elements V, 10c), distorts the outer 
waves 4>. From the upper row in Eq. (0J we find for all 

k 2 ^ £„, 71=1,2,... 



k 2 - Hi, 



-V<p =G- mt V<t). 



(12) 



Here G; n t is the resolvent operator for the internal prob- 
lem. Then, the lower row transforms into 



= 0. 



(13) 



(k 2 - H cxt - G int V 
or, in the position representation, 

(fc 2 + £)#r) + S'(r a+)[0(a) + G^(a, a)0'(a)] = , 
(k 2 + £)0(r) - S(r a + )W{a) G™" '(a, a)0(a)] = 



4 



for the two cases considered. This yields change of the 
BCs of the external part of the exact wave function to 



<P(a) + G^(a,a) 0» = , 
4>'(a)-G ( °l"(a,a) cj>(a) = , 

where the shorthand 

O 2 



Gint («.°) 



3rd; 



(14a) 
(14b) 

(15) 



has been used. Let us note that Eq. I|14a|) is the con- 
ventional boundary condition in the Wigner-Eisenbud R- 
matrix theory of resonance nuclear reactions 0, |28| . 



III. S- AND R FUNCTIONS 

In the external region r > a, where the potential van- 
ishes identically, the wave function has the form <f>(r) — 
const[e" lfcr — S(k)e lkr }. Therefore, the conditions l(T3jl 
obtained above allow us to express the function 



S(k) 



1 



ik<j>{a)/<t>'{a) e _ 2tka 



1 - ik(j)(a)/cj)'(a r 
in terms of the internal Green's functions 



S(k) 



l-ikG™(a,a) 
l + ikG^l(a,a) 



Sf\k) 



l-(i/k)G£l"(a,a) 
l + (i/k)G[f(a,a) 



(16) 

(17a) 
(17b) 



Both expressions are exact and equivalent to each other. 
The only difference is that they are written in different 
complete bases in the full Hilbert space. The merits as 
well as limitations of each of these representations will 
be discussed below. 

We define further the phase shifts <5* N ' D '(fc) due to 
the influence of the internal region and the functions 
i?< N ' D )(fc) ee -2tan<5( N ' D )(fc) by 



RW(k) = 2kG\ N n l(a,a) = 2kj: 



m° (q)M°(q) 



, (18a) 



■«) = ! En 



£'(»-)<(«) ,( 18b ) 



fr 2 - 



The spectral representations of the Green's function are 
used here in the last steps. This shows that the energy 
levels of the corresponding internal problems are real 
poles of the R functions in the complex energy plane. 
The formulas obtained are in close analogy with the rep- 
resentations found in |17| in a different way. They are 
uite similar to those appearing in the i?-matrix theory 
2ST|. We should, however, stress that functions R de- 
fined in our manner differ from the standard Wigner's R 
function by trivial factors like —2k or —2/k which we for 
the sake of convenience include in the definition. 



A certain caution is needed while taking the limit 
r,r' — ► a in Eq. (I18b|) since the derivative of the Green's 
function 



cL D tWo 



9(r> r) X ?(r) X °(r') + 9(r r')x?(r') x°(r) 
^v° n (r)v n (r') 



k 2 



(19) 



has a discontinuity when r = r' . The symbols Xi( r ) and 
%2( r ) stand for solutions of the internal problem with 
Dirichlet BCs only at the points r — or r = a, respec- 
tively. It is readily seen that the mixed second partial 
derivative contains a singular contribution 



X?'(r)x§(r) ~ X°i(r)xl'(r)] 5{r-r') = WS(r 



(20) 



W = — 1 being the Wronskian. This singularity at r = 
r = a must be excluded and the second derivative must 
be understood as 



G i ( n D t ) "(a_, a ) = X ?'( a ) X °'(a). 



(21) 



It immediately follows from this remark that the spectral 
sum ~^2 n Vn'(r) v n ' (r r ) / (k 2 — En ) diverges when r = r' . 
Indeed, convergence of this sum depends on the contri- 
butions of the very high levels, n — > oo. For such an 
excitation, we can neglect the potential U(r) whereupon 
the solution of the internal problem reduces simply to 
v n (r) = ^2/asm.{nir r/a). The contribution 5(r — r') 

to be dropped appears due to the factor — (im/a) 2 
which arises in the numerator after the double differen- 
tiation has been done. It is easy to see, that just the 
contribution of this kind mainly caused numerical prob- 
lems in [13. For r ^ r' , the sum converges, alth ough 
slowly, to the finite limit l|21|) because of oscillations [2!j . 

It would be a mistake to rely upon the spectral rep- 
resentations in (fH^I . and interpret the eigenvalues e n as 
energies of resonance states and the phases £q D ' (fc) = ka 
or Sq (k) = ka+n/2 as smooth phases of the background 
scattering. Indeed, the levels e n depend on what sort of 
boundary conditions have been used. More than that, 
many other forms of BC are equally permissible (see, for 
example, the fundamental review [28j). It can be shown 
that the most general possible BC, 

u°'(a) + I3 int u°(a) = , <f>°'(a) + /WV) = , (22) 

involve two arbitrary parameters (3 (here we include for- 
mally (3 = ±oo for the Dirichlet BCs). Returning to 
Eq. 117JI . we see that, due to the relation Sq (ft) — 
— Sq (fc), the levels £^ D ' are zeros of the function R^(k) 
and, similarly, the levels e n are zeros oiR^(k). In fact, 
the two phases are connected as 

tan (5 (D) (jfe) = - cot 5 {n) (Jfe) = tan[J (N) (Jb) - tt/2] . (23) 

The shift tt/2 just compensates the similar shift of 5^ (k) 
and the total scattering phase S(k) = 5^ (k) + 6^ (k) = 
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S^(k)+S^\k) does not depend on the type of BC used. 
In particular, the positions of the poles of the function 
S(k) in the complex k plane, which are found from the 
equation 

1 - i tan S (D) (k) = = 1 - i tan <5 (x) (fe) , (24) 

are BC independent because of the relation (|23|1 . This 

is in agreement with the fact that both factors sj^'^ are 
entire functions in k plane. 

The factorized form of the residues of poles in i|18|) 
allows one to represent the function S(k) in a different 
but fully equivalent form as 



S(k) = [l-iA 



1 



k 2 -H 



A ) So(k) , (25) 



where the non-Hermitian symmetric operator TL is de- 
fined by 



TL(k)=e - -AA 1 



(26) 



Here £ is a diagonal matrix of the eigenvalues of the 
corresponding internal problems and the column vector 
A(k) of coupling amplitudes, which originate from the 
off-diagonal element V in Eq. J2J|, has the components 



A^(k) = V2ku° n (a) or A^ (k) = ^ ^ {a) (27) 

in the cases of Neumann or Dirichlet BC, respectively. A 
proof of equivalence to the expressions i|17|) immediately 
comes from the following relation between resolvents (see, 
for example, |j|): 

p~h = Gint ~ \ Gint A TTp AT Gint ' (28) 

It follows from Eq. I|25(l that the poles of the S function 
can also be found from the secular equation 



det [z 2 - TL(z)] = 



(29) 



where z is a point in the complex k plane. In the position 
representation this is equivalent to the spectral problem 
with a complex (outgoing-wave) boundary condition 



dr 



r + U{r)j tp(r) = z 2 ip(r) 
tp'(a) — izip{a) = . 



(30a) 
(30b) 



This clearly demonstrates again the independence of the 
choice of BCs for disconnected internal and external mo- 
tions. 

Actually, the form of the operator (|26|l is provocative. 
One is tempted to interpret this operator as an effective 
Hamiltonian whose k-dependent eigenvalues define com- 
plex energies of metastable resonance states formed by 
the potential U(r). However, this interpretation is, gen- 
erally, wrong. To make our points clearer, let us start 



with the following simple remark. Let the potential van- 
ish identically everywhere, U (r) = 0, so that the 5* func- 
tion S(k) = 1. Nevertheless, neither the operator TL nor 
the function 



S(k) = l-iA 7 



1 



k 2 -H 



A 



(31) 



are trivial in themselves. For example, in the case of 
Neumann BCs matrix elements of TL are 



H W = JL 

or 



7T 2 <5„ 



ika{-l) {m+n) 



(32) 



Any truncated finite-size N x N matrix obtained from 
Eq. I|32(l gives, when substituted in Eq. H29J) . N pairs 
of complex roots. However, the complex poles of the 
resolvent [z 2 — 7itrunc(^)] _1 of the truncated matrix have 
nothing to do with the poles of the genuine S function. 
The truncation procedure is not stable when N — > oo. 
The difficulties become even worse in the case of Dirichlet 
BCs when the imaginary part also grows, 



n 



(D) _ 



1 r 
~7p- 



ka 



(-1) 



(m+n) 



(33) 



In fact, all poles found go to infinity in the limit N — > 
oo. Indeed, the true poles, because of the identity 
det[z 2 - H(z)] = det(z 2 - H int )[l + %R(z)], must sat- 
isfy the equation 



1 



-R(z) 



0. 



(34) 



which is equivalent to l|24|) . For example, in the case of 

Neumann BC this equation runs as 1 + iz G[^ t ' (a, a) = 0. 
In particular, when the potential U(r) vanishes iden- 
tically, the Green's function in the z plane is equal 
to G^2i a ,a) = -(l/z) tanza, Eq. looks as 1 - 

i tan za — and, therefore, has no roots in any finite 
domain of this plane. This implies, in turn, that Eq. I|30j) 
has no non-trivial solutions. In fact, the function H31(l is 
simply equal to Sf Tee (k) — e 2lka in this case and compen- 
sates exactly the phase shift due to the fictitious reflec- 
tion at the separation point a. 

The remark above is of quite a general nature. The ra- 
dius of separation can be chosen arbitrarily. For the sake 
of simplicity, we suggest only that a is larger than the 
finite radius of the potential. The actual wave function 
satisfies at this point the conditions 



(f)(a) = u(a) , 4>'{a) = u'(a) . 



(35) 



Of course, none of these quantities is known before the 
problem has been solved. Any boundary condition used 
above generates a complete basis in the Hilbert space, 
in which the actual wave functions can be expanded. 
All such bases are formally equivalent. But this does 
not mean that all of them are equally adequate from the 
physical point of view. In particular, the more the ba- 
sis eigenvectors and their derivatives differ at the point a 



6 



from the real values i|35|l the more slowly the correspond- 
ing expansion converges near this point. 

A forced BC creates a false reflection at this point, 
which is described by the factor So(k), whose phase 
should be fully compensated by the similar part of the 
total phase of the function S(k). Both factors separately 
depend on a although the complete function S(k) does 
not. The main role in this compensation is played by the 
matrix elements TL m n with large m, n when the influence 
of the potential U becomes negligible, and we return to 
the situation described in the previous paragraph. 

However, the choice of the separation radius and BC 
influences the positions of the poles (as well as the 
residues) of the function R and, consequently, the ex- 
plicit form of the matrix Tl mn . This influence is stronger 
the less adequate the choice made of BC and a. In gen- 
eral, the parameters of the R function can carry rather 
poor information about the actual complex poles. For 
this reason, diagonalization of a truncated matrix Tl mn 
(which is necessary in any numerical computation) can 
lead in the case of a poor choice to strong disagreement 
with the characteristics of the actual poles of the function 
S. The explicit dependence of the matrix elements TC m n 
on the wave number k causes additional problems (see 
next section). To extract the physically relevant effec- 
tive Hamiltonian from the formal operator 7i, additional 
physical considerations must be engaged. For example, 
one may expect from the physical point of view that the 
most relevant choice of the separation point a would be 
a distance matching an outer potential barrier which is 
strong enough to make immediate reflection at this point 
quite probable. 

To explore in more detail the questions briefly dis- 
cussed above, we will apply in the next section the formal 
technique sketched here to the problem of scattering by 
a finite periodic set of 8 barriers, which can be solved 
exactly. 



IV. OPEN KRONIG-PENNEY MODEL 



tion, the equivalent integral equation 

/>oo 

X (r) = sin(fcr)+ / dr'G° + {r,r')U(r') X (r') 
Jo 

L 

= sm(kr)+Y,KiG° + (r,l) X (l)- (37) 
i=i 

Here the symbol G+(r, r') — — (1/fc) g°(r, r') with 
g°(r, r 1 ) = 6(r' - r) sm(kr)e lkr ' + 9(r - r 1 ) sm(fcr')e ifer 

stands for the Green's function of the free radial motion 
which has outgoing-wave asymptotic. From the second 
line in 1|37|1 we find immediately S(k) — 1 — iT(k) where 
the scattering amplitude is given by 

T ( fc ) = ^E^ sin ( fcZ )^)-^ T « x . (38) 
i=i 

In the second equality we have used matrix notation, with 
s and x being L-dimcnsional vectors with the compo- 
nents si = si(k) = sin(fcZ) and \i = x(0 (£=1,2,... ,L), 
respectively, when k = diag{K l5 k 2 , 

According to Eq. I|37|) , the L-dimensional vector \ sat- 
isfies the equation \ — s + G a n\, where G° = —(l/k)g° 
and g° is a symmetric non-Hermitian matrix with the 
matrix elements 



9°w 



sin(fcZ) e m ' if I < V . 
sin(M') e lkl if I > V . 



(39) 



Thus, we obtain finally 



S{k) = e 2l5 « = 1 - z S T S(fc) S = 1 ^ K [ k \ , (40a) 



1 + ±K(k) 



K{k) = -2tan6(fc) = s T G(k)s . 



(40b) 



where we have used the factorized form 3m g° — ss T of 
the anti-Hcrmitian part of the matrix g° to pass from the 
first to the second equality. The LxL matrix propagators 



A. Exact solution 

We will consider below s-wave scattering by a period- 
ically disposed chain of L 5 barriers, 



0-(r)=J>*(r-O 



(36) 



i=i 



To ensure formation of long-lived resonance states, at 
least some of the strength constants ki should be posi- 
tive. The distance r is measured in units of the size of the 
well formed by two neighboring barriers. Below, we de- 
rive the effective Hamiltonian starting directly from the 
Schrodinger equation. 

Because of the local character of the barriers, it is most 
convenient to start from, instead of the Schrodinger equa- 



9{k) 



k\ + g° 



and G{k) 



kX + 3?e g° 



(41) 



are connected with one another by a relation similar to 
B28J1 . The diagonal matrix A = ft -1 characterizes the 
penetrabilities of the barriers [3(J. The poles of the S 
function are defined by the equation 



det[2A + g°(z)] = = 1 + | K(z) . 



(42) 



All the matrix elements Q39JI are entire functions in the 
complex z plane. The same is valid for the determinant 
in H42|) . Therefore, one can show that this equation has 
an infinite number of isolated complex roots. Further, its 
roots because of the relation [^/(z)]* — ~9u'(~ z *) come 
in pairs z n and Z- n = — z*, n=l,2,..., symmetrically with 
respect to the imaginary axis, or lie on the latter, z q = 
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iy q . All poles of the first type are situated in the lower 
part of the complex plane. Poles on the positive half of 
the imaginary axis, y q > 0, correspond to bound states 
and can appear only if some number of the constants k; 
are negative. Those which are situated on the negative 
part, y q < 0, correspond to the so-called virtual levels. 
The total number of purely imaginary poles is finite for 
any finite L. 

After all poles have been found, the function S(k) can 
be presented in the form 



s(k) = n 



(fc + z n )(k - z*) 
\ Zn){k + Z*) 



n 



fc- 



Wq 



(43) 



where we took into account that S(k = oo) = 1, since S 
barriers become transparent for a particle with asymptot- 
ically large energy. These expressions are in agreement 
with the general theory |3ll l32]| , which is valid for any 
potential with a finite radius. 

Each factor in (|43[) is singly unitary. In particular, for 
a pair z± n of conjugate roots we have 



S„(k) 



_ (fc + z n )(k - z*) 
(k-z n )(k + z*) 



2ifc3m Zm 



fc 2 - 



(44) 



Since Smz„ < 0, the phase S n (k) increases when the 
energy E = fc 2 grows and passes the value tt/2 at the 



point E 



|2 — 



E n . The typical energy interval AE 



of the main gain of the phase is estimated by the quantity 
2fc|3mz n |. If this interval is small enough and the latter 
quantity does not vary appreciably within it, the total 
gain is close to it and the factor S n receives the standard 
Breit-Wigner resonance form 



E — E.„ 



-T 

2 n 



E-£* 



E — E n + iT n E — £ n 



with the energy E n = $ie£ n and width r„ 
of the resonance defined as follows 

E n = \z n \ 2 , r„ = A\1Rez n Qmz n 



(45) 



-23m £ n 



(46) 



The fc-dependence is neglected in the definition of the 
width and substitution fc = \f~E~ « |9?e z n \ has been made. 
Such a substitution is well justified when the scattering 
energy is large enough, but becomes improper near E — 
0. Due to the "threshold" \/E dependence of the widths, 
some specific behavior takes place when in the proximity 
of this point there exists a number of bound and decaying 
states [H Ulil. 

We conclude that the scattering amplitude T(fc) = 
i[S(k) — 1] is a meromorphic function in the z plane and 
can therefore also be written down as an infinite sum of 
the pole contributions 



Qn 



Qq 



(47) 



where the residues g n:q can easily be found from Eq. I|43|l . 

Another representation of the scattering amplitude 
sometimes considered in the literature ensues from diago- 
nalization of the (finite in our case) matrix Q (fc) at a given 
fixed real fc. Since this matrix is complex symmetric, such 
a diagonalization is performed by a complex orthogonal 
transformation defined together with the complex eigen- 
value matrix A(fc) = diag{Ai(fc), A2(fc), A^(fc)} from 
the equation 



[fcA + g°(fc)]*(fc) = #(fc)A(fc). 



(48) 



The scattering amplitude reduces then to a finite sum 



T(fc) = 2E' 



i=i 



s T (fc)^W(fc)] s 
Mk) 



(49) 



where vE'^(fc) is the Ith eigenvector. However, such a 
representation is of limited practical use inasmuch as the 
terms of the sum are, generally, extremely complicated. 
There is no simple interpretation of a given term and all 
of them can give comparable contributions. The com- 
patibility of the two representations l|47() and l|49(l is very 
indirect. A particular term of the second one cannot 
be uniquely continued from the real axis at an arbitrary 
point z in the complex plane because both ^^ l \z) and 
A;(z) are multivalued functions in this plane. The same 
is, of course, also true for the opposite direction. There- 
fore, even if a root zj of the equation A;(z) = is found 
and near this point the pole approximation 



T\{z) 



2[,s T (z^ W (4)] 2 /A^ 



(50) 



is valid, this expression cannot, generally speaking, be 
continued on the real axis by simply substituting z — » fc, 
since the power expansion Ai(z) — A;(zj) + AJ(z')(z — 
z l j) + ... has a finite radius of convergence. Further, for 
any fixed I a set of roots exists depending on the branch 
of the function A;(z) considered. There is no one-to-one 
correspondence between the set of roots zj and the man- 
ifold of actual poles z„ i9 that appears in the expansion 
l|47p. Many of the roots zj are false and their contribu- 
tions must finally cancel out. Therefore, diagonalization 
of the propagator with non-trivial energy dependence is 
not as a rule useful; rather it can lead to misleading con- 
clusions. 

We compare below the closed expressions Eq. (|40|) 
found with the formalism described in the previous sec- 
tion. In correspondence with our remark at the end of 
Sec. Ill, we fix the separation point a by superposing it 
on the position of the outer barrier, a = L. The appear- 
ance of an additional S function changes the correspond- 
ing boundary conditions. In particular, if the <5-terms are 
ascribed to the internal region, Neumann BC J7|) for the 
internal problem is replaced by 



K L u°{a) + u°\a) = = u°{a) + X L u ' '(a) , 



(51) 
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whereas attributing such terms to the external domain 
yields the following change of the BC of the external 
problem: 



K^{a) - <P°'(a) = = <f>°(a) - A L ^'(a) . 



(52) 



We stress that neither of these conditions coincides 
with the boundary conditions 



X(o-) = X{a+) = x(a) , 
X'{a+) - x'(a-) = KlX(o) 



(53) 



satisfied by the exact wave function x( r )- 

In the case (|52ll . the immediate reflection at the point 
a is described instead of by 



S (k) = S' a 



(N)' 



(Recall that S^\k) = er 2lka = e~ 2lkL .) The additional 
factor is due to the influence of the barrier outside the 
radius of separation a. This factor has a pole k = —ik l on 
the negative part of the imaginary axis, which formally 
corresponds to a virtual Wigner level. Actually, such a 
pole does not exist in the exact solution and disappears 
due to cancellation with the contribution of the internal 
region [see Eq. (|95|l below]. 



with the levels £n being the eigenvalues of the internal 
problem with the BC l|5ip. For the sake of simplicity, we 
use the same superscript N as before. 

In spite of the seemingly similar general structure of 
the expressions (|40|) , <|41[) on one hand side and l|55[l 1571) 
on the other, they are, in essence, quite different. The 
most important distinction shows itself in the dimension 
of the vectors and matrices which is finite and coincides 
with the number of barriers in the first case and infinite 
in the second. In addition, the factor So(k) of the im- 
mediate reflection does not appear explicitly in Eq. I|4(J|I . 
We will analyze below a couple of simple special cases 
before drawing general conclusions. 



C. One 5 barrier 

In this case we find immediately from Eqs. (|40|l (a = 
L = 1 and we drop the subscript L in the strength of the 
barrier) 



S(k) 
K(k) 



sin ke %k + A k 
sin k e lk + A k ' 

2 sin 2 k 
sin k cos k + Afc ' 



(60a) 
(60b) 



whereas 



B. Internal problem with Neumann BC 

The exact S function in the case of the BC (|51|l reads 
as S(k) = S(k)S^\k), where 



with 



R(k) = A 7 



\2 V^°° U n ( a ) U n ( a ) 



and 



f(k) = A 1 



1 



k 2 - 



A. 



(55) 



A = 2kG\l>(a,a) (56) 



(57) 



The coupling amplitudes are equal to 

A n {k) = V2ku° n (a) = -\ L V2ku° n '(a) . (58) 

Finally, the matrix elements of the operator 7i' N ' [see 
Eq. appear as 

Hmn ~ e ^m &mn — ik U° m {a) U° n (a) 

= e^5 mn -z\lkv?J(a)^(a), (59) 



sin k + Afc e 



ik 



sin k + Afce %k ' 

„ , , x „ , k sin k 

R(k) = -2A— 



sin fc + A k cos fc 



(61a) 
(61b) 



In the Eqs. (|61|l the phase 5^(k) = —k of the immediate 
reflection at the point a is extracted. The two functions 
K{k) and R(k) are related to each other by 



K(k) 



2 tan fc + R{k) 
1 - i tan fc R(k) ' 



The positions of the poles of the function S coincide 
with the complex roots of the equation 



+ 2i\z~- 1 = 0. 



(62) 



In the right half of the z plane they can be searched for 
in the form z n = nn + £„, where any is restricted to 
the strip |3?eCn| < tt/2 and satisfies the equation 



e„ + -m[l-2iA(mr + e„)] =0. 



(63) 



There also exists the trivial root z — but this root 
fully cancels out. Even for not very large n one can omit 
C n in the logarithmic term. This yields the approximate 
solution 



ln(l - 2i\nir) 



(64) 
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which is valid with good accuracy for almost all poles, 
being asymptotically exact when n — > oo. 

However, by no means all of the complex roots z n cor- 
respond to resonances. For a pole to correspond to a 
long-lived resonance state, the following two additional 
conditions also must be satisfied: i) 5Rez„ should be a 
real pole of the function K (k) or at least should be close 
to such a pole; ii) the corresponding residue should be 
small enough for the scattering phase S(k) to increase 
near this point by ir in an interval Afc <C 1- The validity 
of these conditions depends on the value of the parameter 
A. Indeed, the poles of K(k) are found from the equation 



sin k cos k + Xk = . 



(65) 



There exists only a finite number of real roots ±kj, 
j=\,2,...,j max of this equation. They satisfy the require- 
ment 2Xkj < 1 or, equivalently, kj < k/2. (The trivial 
solution fc=0 should be dismissed.) All other roots lie in 
the complex plane and the phase does not reach its maxi- 
mal value although can change rather fast if a pole of K is 
still close to the real axis. It is clear that they move away 
from this axis when \z\ grows. In particular, no real solu- 
tions exist if A > Ao ~ 0.2 (k < n m 5). The scattering 
phase S(k) < ir/2 in this case and smoothly depends on 
k. The barrier is too weak to form a long-lived resonance 
state. When «| sinfc/fc| < 1 the scattering phase can be 

calculated as S(k) w —2k sin 2 k/k-\ in the framework 

of perturbation theory. 

There are two kinds of real roots of Eq. I|65|) when A 
is appreciably less than Ao- The most interesting case 
Xkj <C 1, although j ^> 1, can easily be considered ana- 
lytically. The first set consists of the roots k n which are 
close to the points nir: k n rs (1 — A + X 2 )nir, (nir < 1/A). 
Near such a pole the K function manifests typical reso- 
nance behavior 



K{k) 



2{Xmif 
k h n 



(66) 



Strictly speaking, a similar contribution of the symmetric 
root k- n = — k n , which corresponds to the same energy 
£n = fc 2 , should be added, so we arrive near the nth 
resonance at 



K n (E) 



(2X) 2 (mrf 
E — £„ 



(67) 



Each neighboring pair of resonance roots is separated by 
a root of the second set: k m w [1 + (Am7r) _1 ](m + 
in a vicinity of which 



K(k) 



(68) 



The residues are large in this case and such terms con- 
tribute into the smooth background part of the total 
scattering phase. Indeed, the (dimensionless) time de- 
lay T(k)—d5(k) / dk=r w /\tq\, which measures the Wigner 




FIG. 1: One S barrier model for A=l/35. (a) The time delay 
T(k)=dS(k) /dk, (b) K function and (c) ii function. Only 
a finite number n rcs ~ [l/27rA] (= 5) of poles located near 
the points Trn, n=l,...,n ICS , correspond to long-lived resonance 
states; the others refer to the smooth phase of reflection at 
the point a(=l). 

time delay T w =2d5/dE in units of the (negative) time 
delay To=—2a/(dE/dk) because of the immediate reflec- 
tion at the point a(=l), is large near the points k n but 
small when kK,k m . Poles of the two different kinds alter- 
nate and resonant and smooth contributions are mixed. 
Figure 1 clearly demonstrates all the features described 
above. 

Let us on the other hand consider the properties of the 
function R(k) [see Eq. I|61|)]. Its poles are the roots of 
the equation 



sin k + X k cos k = 



(69) 



instead of Eq. I|65|) . Contrary to the latter, all roots ±k n 
of (|69|l are real. Those of them, that satisfy the inequal- 
ity Xk n <C 1 and, therefore, have numbers n < n max , are 
equal to fc n w (1 — A)n7r again and correspond to reso- 
nances. The difference from the similar roots of Eq. I|65|) 
appears only at higher order corrections in the parame- 
ter A. For larger n, when Xk n 3> 1, an infinite number 
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of roots k n « ir(n + |)[1 + A 1 (nir) 2 ] exists, giving 
a smooth contribution, which combines with the phase 

S^\k) — —k and almost compensates it at large k [see 
Fig. 1(c)]. Indeed, the S barrier becomes almost trans- 
parent when the collision energy is large and the total 
scattering phase can be calculated by applying perturba- 
tion theory. 

The meromorphic character of the function R(k) al- 
lows us to represent this function in the form of the pole 
expansion 



*(*) = E 7 



(Xkj, 



n— — oo 



- 2fc E T 



+ A + (Afc„) 2 k - k n 
2 (Afc„) 2 



\ 1 + A + (X~k n ) 2 k2 - 



(70) 



The symmetry connection fe_„ = —k n has been taken 
into account in the second equality. 

Returning to the poles of S(k), we can expand the 
expression (|64|) in the resonant region 2Xmr <C 1 with 
respect to this parameter to calculate the poles 



z n m (1 — A + A 2 ) mr — i(Xnir) 2 
and the complex energies of the resonance states 
En « (1 - 2A + 3A 2 )(n^) 2 - |[4A 2 (n7r) 3 ] . 



(71) 



(72) 



Note that the widths T n = 4A 2 (rt7r) 3 of the resonances 
appear only in the second order in the penetrability pa- 
rameter A when the shifts of their energies are of the first 
order of magnitude. The resonances are well isolated 
since the ratio T n /(E n+ i — E n ) = 2n(Xn) 2 -C 1. Finally, 
the remote poles with n ^> l/27rA are given by 



z n w (n — j)tt — | ln(2An7r) 2 



(73) 



Now we compare our findings with the results l|55[) - 
()57|l of the general formalism. The normalized solutions 
of the internal problem with BC l|51|) are readily found 
to be 



1 + A cos 2 k r 



■ sin(fc n r) 



(74) 



where k n are the roots of Eq. (|69|l which follows directly 
from the BC (|51|) this time (we omit the tilde and note 
that only positive roots k n , n.=l,2,..., are to be kept). 
This fact is quite satisfactory and demonstrates the phys- 
ical relevance of the choice of BC made. Indeed, for such 
a BC u n (a) ~ A in accordance with the exact boundary 
condition JHJ). The Dirichlet BC u° n {a) = would be 
deficient in this sense. With Eq. I|74|l taken into account, 
we find from Eq. (I5C 



i?(fc) = 2fc]T T 



2cos 2 fc„ (Afc„) 2 



+ A cos 2 k n k 2 — e„ 



(75) 




FIG. 2: Widths (in log scale) versus energies for A=l/35. 
Exact poles [from Eq. 16311 . □] and approximation [Eq. 1721 . 
0], compared to the "diagonal" approximation TC nn (k = k n ) 
with Neumann [from Eq. I76L x] or Dirichlet BC [from Eq. 
I|97|l . +]. Inset shows the asymptotic region of remote poles, 
where the diagonal approximation breaks down. 



Equivalence to Eq. (|70Jl is seen from the relation 
cos 2 k n = [1 + (Afc„) 2 ] -1 , which follows directly from the 
secular equation (|oT)|) . 

Matrix elements of the operator Ti are in our case 



"17 (*0 _ px 

1 l"mr> ^<n u T7 



-i2k- 



Xk r , 



Xk r , 



yl + A + (Afc m ) 2 v /l + A+(Afc„) 2 

(76) 

In the resonance domain, Xk n < 1, the diagonal matrix 

element Tinn(k — k n ) approximates well the complex en- 
ergy £ n of the resonance which lies near the scattering 
energy E ss k 2 . In this region off-diagonal elements influ- 
ence only corrections of higher orders. At the same time, 
the diagonal approximation becomes insufficient for the 
remote poles, breaking down to reproduce exact asymp- 
totic nlnn behavior of the imaginary parts. Figure 2 
illustrates the consideration. 

We therefore conclude that the extraction of the phase 
<5q D ' (fc) of the immediate reflection and utilization of the 
R— function is adequate only in the resonance region 
where such a reflection is probable. Beyond this region 
the internal and immediate reflection phase shifts almost 
compensate each other so that the function K (k) proves 
to be a more relevant tool. 



D. Resonance domain 

Below, we restrict our attention to the resonance re- 
gion well above the energy E — k 2 — 0. To simplify the 
investigation further, we suggest assuming the strength 
constants to be the same A; — X\, I — 2,3, ...,L— 1, for 
all inner barriers. In accordance with the results of the 
simple consideration presented at the end of Sec. IV C, 
segregation of the smooth phase of the immediate re- 
flection meets expectations based on physical intuition. 
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To perform this segregation, we write first the scattering 
amplitude T{k) — s T Q{k)s as a block product 



T(k) 



Q T 

(A L fc - s L e ifca )s T 



1 



X^I + g 



(77) 



where the symbols Q and etc., stand for [L — 1) x 
(L — 1) submatrices and (L — l)-climcnsional vectors, re- 
spectively, and g° is the similar upper left block of g°. 
To pass from the first to the second line, the relations 



s T Qs = - (X L ke 



-ika 



and 



s T F = ;f t s 



_ e ika s T 



1 



\ikl + g° 

have been used which follow, together with the expression 



Gl = 2 [\ L kc 



-ika 



e lka s T 



1 



Xikl + g 



■si 1 e~ lka . 



from the equation (kX + g°)Q(k) — 2. We notice now 
that the function 



complex poles of the S function well separated from all 
the other poles. 

The corresponding R function is readily found from 
Eq. O to be 



R(k) = - 



1/X h k + cot(<5oo + ka) 



(81) 



In the limit A L — ► (closed interior) we have R(k) = 
and S(k) = 1, so that only the immediate reflection at 
the point a survives. Obviously, the spectrum of the poles 
of the function R(k) is determined by the equation 



sin (Soo + ka) + X L k cos (Soc + ka) = 



(82) 



which should be compared with Ea.! f69p . It is easy to 
check that the spectrum exactly coincides with that of 
the wave numbers of the internal problem for the poten- 
tial Uoo(r) = Ki Y^=i S(r - with the BC (£>TJ. This 
potential is perfectly transparent on the separation radius 
a = L where the latter condition is fixed. This directly 
follows from the expression 

tt°(r) = sin(fcr) - \ T^k) e lkr = § [e~ lkr - S^k) e ikr ] 

(83) 

for the wave function in the region (L — 1) < r < a = L. 



Xikl + g° 



(78) 



coincides with the amplitude of scattering on the poten- 
tial with the last barrier being removed: k l = or, 
in other words, A L = oo. Finally, we obtain after simple 
transformations S(k) = S(k) e~ 2lkL , where 



S(k) 



-2ikL 



(l + 2iX L k)S oc (k) 

(l-2^A L fc)e- 2 ^-5 0O (fc) 
sin(<5 00 + ka) + X h ke i( - s ^ +ka '> 
sin(£oo + ka) + X h ke- i( - s °°+ kaS > 



(79) 



The poles of the function S(k) [as well as of S(k)} in 
the complex k plane are given by zeros of the denomi- 
nator. It is convenient to introduce in parallel with this 
function a sequence of functions S®(k) = S^(k)e 2ilk 
which describe the scattering on chains of I barriers 
with the penetrability constant Ai and the BC fixed 
at the radius a = I. In particular, S^(k) = S(k), 
S( L -V(k) = S 00 {k)e 2li - L ^ k , and S^(k) = 1. These 
functions are related to each other by a recursion 



S {l \k) 



l-(l + 2iAifc) e 2ik SV-V(k) 
\-2iX 1 k-e 2ik S^- 1 ){k) 



(80) 



Using as before the ansatz z m = mix + £ m in the res- 
onance region Ai, L |z| Ai, L |m|7T <C 1, we arrive at an 
algebraic equation p( L )(e 2 ^ m ) = 0, with V(w) being a 
polynomial of the Lth power with respect to the argu- 
ment w = e 2 ^" 1 . This equation gives a bunch of L close 



E. Two 5 barriers; resonance trapping 

Now we will use the formulaes just found to analyze 
as an illustrative example resonances in the double-well 
^-potential, L — 2. The equation for the poles of the S 
function Ij79(l looks in this case like 



(1 + 2iXiz) e 4iz - 2(1 - 2iX 2 z) e 2lz 
+ (1 - 2iAiz)(l - 2iX 2 z) = 0, 



(84) 



[cf. Eq. H62(l ], Again, the ansatz z m = mm + Cm with 
l^eCml < 7r/2 is valid. Supposing also the integer num- 
ber m to be appreciably large, we arrive for each fixed m 
at a couple of closed solutions [compare with (|64J) ] 



ln<M 



. (2Ai + X 2 )mir 
1 + 2zAim7r 



(85) 



1±4 1- 



4AiA 2 



(2Ai 



(1 + 2iAim7r) 



It is worth noting that within the approximation adopted 
the sum of the imaginary parts 3m (£+ + £~) = 
-| ln\Z+Z~\ = -| In y/l + 4Xl(mir) 2 , where Z± stand 
for the arguments of the logarithm in Eq. 1)8 5[l. does 
not depend on the penetrability constant Ai of the in- 
terior barrier. Indeed, simple transformations show that 
Z+Z- = [(1 - 2iAim7r)/(l + 2iX 1 mn)} (1 - 2iX 2 mn). 

We fix now the number m 3> 1 within the resonance 
domain, Ai 2|Ti7r| 1, and consider the doublet of poles 
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close to this point. Keeping the terms of the two first 
orders of magnitude we find 



1 



Cm ~ ( 2Al + A2± V 4A i + A 2 1 " l7i 



-An 1 ± 



V4A 2 + A 



(mir) 2 . (86) 



In the resonance region, the splitting |Az„ 



I AG 



yl 4A| + A| within a given doublet is much smaller than 
the distance between adjacent doublets which is ~ ir. 
Equation (|81|l reads in the case L = 2 

2A 2 fc (sin 2 k + A x fc sin 2k) 
sin 2 fc+ i(2Ai + A 2 )fcsin2fc + AiA 2 fc 2 cos2fc ' 

(87) 

For the mth resonance doublet k m — mir + 8k m , and the 
small shift 5k m satisfies the quadratic equation 

(Sk m ) 2 + (2Ai + A 2 )(m7r) Sk m + AiA 2 (m^) 2 = , (88) 

giving immediately 5k^ = SfteCm, with from (JSSJl. 

This convinces us that each resonance doublet with a 
good accuracy can be considered independently of the 
other poles. Near a given k m = mir within the resonance 
region the function R(mir + Sk) = R m {k) is decomposed 
into a sum of two partial fractions 
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1 

in 
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0.01 
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X 
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x X 




E m -(mn) 

FIG. 3: Widths (in log scale) versus energies for a doublet of 
resonances with m=10 as (a) A1/A2 is varied from 0.1 to 10 
(marked with •) at fixed Ai=0.001; (b) A1/A2 is varied from 
1 to 10 (marked with •) at fixed A 2 =0.001. 



R m {k) = 



17, 



In 



(89) 



the residues being found to be equal to 7^ = — 23m 

To pass from the k to the energy plane, the term 
R-m{k) must be added, which yields 



Rm(E) = 



r+ 



E - el 



E - el 



- A 1 



1 



E-e 



(N) 



An 



(90) 



The two-dimensional vector A m =(\/r m , y/T m ) and ma- 
trix e™'=diag (e+ , e m ) of the internal levels e„=(fc^) 2 
of the doublet have been introduced on the last step. 

The corresponding representation of the scattering am- 
plitude near the scattering energy E w fc 2 ^ is as follows: 



f m {E)=A\ 



E 



\ / 7 



■ A r , 



(91) 



In these formulas = 2yE^^ w 27rrn7^, i.e., we have 
neglected the change of the scattering energy E within 
the doublet considered. The difference between fc^ con- 
tributes only at higher orders. This is contrary to the 
case of widths from different doublets when k m differ 
from each other already in the zero order in Ai j2 . In such 
an approximation, both functions R m (E) and T m (E) be- 
come meromorphic in the complex energy plane. 



By the very construction, the complex eigenvalues of 
the two-dimensional energy-independent symmetric ma- 
trix 



n. 



(N) 

off 



,(N) 



AA 1 



(92) 



coincide with the complex energies of the two resonances 
of the doublet considered. Henceforth, we drop the in- 
dex m of the doublet to avoid too bulky notation. The 
energies e^ are the levels of the internal problem for 
the potential U~£\r) = Ki5{r - a) with the BC l(5T|) 
whereas the amplitudes are given by Eq. H58|) . with the 
wave number k substituted by 77171". Just the energy- 
independent matrix l|92|l is naturally interpreted as an 
effective non-Hermitian Hamiltonian. The notion of the 
effective Hamiltonian is valid, however, only within a 
fixed doublet. Similarly, in the case L = 1 the "effective 
Hamiltonian" coincides with the corresponding diagonal 
matrix element of the formal operator 7i constructed in 
Sec. III. 

As we have already mentioned above, the total width 
r = r + +r_ « 4A 2 (TO7r) 3 of the doublet remains constant 
as long as the openness of the system is fixed. At the 
same time, the individual widths depend also on the ratio 
Ai/A 2 . A remarkable redistribution of the total width 
between members of the doublet takes place when this 
parameter changes. When the system is almost closed, 
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A 2 < Ai, 



E+ ps (1 - 4A a - A 2 ) (mTr) 2 , 
ps (1- A 2 )(m7r) 2 , 



:2A^(tott) 3 . (93) 



Both levels have similar widths in this case. However, in 
the opposite limit A 2 3> Ai the state that exists in the 
outer well appropriates almost the whole total width 



E + ps (1 - 2A 2 - 2Ai)(m7r) 2 . 

w (1- 2Ai)(m7r) 2 . 



r. 



4(A 2 -A 2 )(m7r) 3 , 
4A 2 (m7r) 3 . (94) 



Simultaneously, the energy E + of the broader resonance 
gets strongly displaced due to the coupling to the energy 
continuum. The other resonance turns out to be trapped 
in the inner well. Figure 3 illustrates the behavior just 
described. For an arbitrary ratio Ai/A 2 , the complex en- 
ergies of the two resonances of a doublet are separated by 
a distance which is large on the scale of the total width 
|35j . For this reason, the off-diagonal matrix elements 
of the anti-Hermitian part of the effective Hamiltonian 
(|92|l give only corrections of higher order and can there- 
fore be neglected. In spite of the nontrivial behavior of 
the complex levels just described the resonances may still 
be considered separately of each other within the given 
accuracy. 

Extension to the general case of an arbitrary number 
L of the barriers is now straightforward. The effective 
Hamiltonian appears as a L x L block of the operator 
7i, which describes a fixed bunch of L close resonance 
states. It is important that energy dependence can be 
fully neglected within such a bunch. 



Comparison with the closed solution (|55|l and (|79(l shows 
that 



1 (d)" 

cot (6oo + ka) = -^G\ n t ( a > a ) 



1 ^<VKV) 



k U2 _ ' 

(96) 

which implies that the spectrum of the internal problem 
is given by the equation sin (5oo + ka) = instead of 

Eq. (|82(l . Unlike the eigenvalues Sn \ the levels do 
not depend on A L . 

Similar to Eqs. lj55j) - lj57|) . we have S(k) = S(k)S (k), 
with S from Eq. O, and S(k) = l-iA T [k 2 -H^^A, 
where 



A, 



n mn s m o mn 1 _ iX , 



W m'(°) V n( a ) 



e iD) S 



1 - - i - - 
2^ ^ A m A n — — A m A n (97) 



and 



A, 



2k 



A 2 /c 2 



v° n \a) 



(98) 



This should be compared with Eq. and lj59"j) . Not 

only widths but also Hermitian shifts result from the 
coupling to the external region this time. Indeed, in the 
resonance region X h k <C 1 the phases that come from the 
first two factors in Eq. i|95b|) are equal approximately to 



Gtl (a, a) 



8(k) ps -A 2 fc- 



A L l_r int 



(a, a) 



F. Internal problem with Dirichlet BC 

The internal problem with the BC l|51(l does not corre- 
spond to a finite motion of a quantum particle even when 
the coupling matrix elements are neglected. Indeed, the 
internal solution is sensitive to the penetrability of the 
outer barrier. Meanwhile, in numerous applications the 
intuitively most attractive convention, which goes back 
to the textbook 0, is adopted, using as an internal basis 
the states of a closed counterpart of the system under 
consideration. Such a counterpart can hardly be defined 
uniquely but, in our case, it is natural to fix it by choos- 
ing Dirichlet BC at the point r — a. This corresponds to 
the internal problem with an infinitely hard wall put at 
the point r = a. The formalism of Sec. II gives then 



S(k) 



1 + + (a, a) 



S, 



(N) 



1 — i\ L k 



The two contributions perfectly compensate each other at 

the points fc 2 = . Within a bunch of resonance levels 

the eigenvalues of the Hermitian part of coincide 
in the main approximation with the corresponding levels 

■ In general, the connection S(k) = 5(k)+S\ L (k) holds 
for arbitrary k. 

Let us first return to the case L — 2, for a moment. The 
amplitudes (|98|l for the mth doublet are easily calculated 
to be in the main approximation A± ps ±V2A 2 (m7r) 3 / 2 , 
when the corresponding internal levels are e^ ' ps (1 — 

4Ai) (m7r) 2 and = (nm) 2 . Therefore, the diagonal 
elements of the effective Hamiltonian of the doublet are 
equal to 



4 0) 



(1 - 4Ai - A 2 ) (mTr) 2 - i\ 2 2 (mn) 3 , 
(1 - A 2 ) (mTr) 2 - iArJ(mTr) 3 , 



(99) 



(95a) which exactly coincide with Eq. I|93l) . Off-diagonal ele- 
ments of can be neglected if A 2 <^ Ai. These ele- 
ments become important, however, when the coupling to 

(95b) the continuum becomes strong, A 2 > Ai. The two reso- 
nances (|99H interfere in this case. Diagonalization of the 



14 



2 x 2-matrix Ti^a gives now 



£± 



1 - 2Ai - A 2 t 



4A 2 



1 ± 



(mir) 



(100) 



In the limit A 2 >• Ai this reduces to <|M|) . 

For an arbitrary number L of barriers, the mth bunch 
of L close resonances is described by the L x L-submatrix 
7igg of the infinite matrix l|97l) . with the wave number 
fc being substituted by k m ss run. The absolute val- 
ues of the coupling amplitudes within such a resonance 
bunch are estimated as \A n \ ss 2A L (m7r) 3 / 2 / \[L. The to- 
tal collective width if' = Ya=\M = A 2 ~ 4A 2 (m7r) 3 
of the bunch, which is determined by the trace of the anti- 
Hermitian part of "H^a , characterizes the openness of the 
system and does not depend on the penetrability con- 
stant Ai of the internal barriers. The same is valid in the 
main approximation regarding the collective real energy 
displacement Sei = 2A L (m7r) 2 w vi" 1 * 1 /2\ L rmr. When 
this displacement is small in comparison with the total 



energy spread of the internal levels, 5e c <C A^ m ^(e 
the bunch consists of L independent narrow resonances 



c(m) 



- 8e c "^ 
L 



i 1 
2L 



p(m) _ 



(101) 



Note that the Hermitian shift does not influence the level 
spacings but changes only the position of the bunch. 

Under the opposite condition, Sec ^> A' 1 ™) [s^ D })_, the 
doorway basis in which the interaction matrix AA T is 
diagonal becomes more adequate j2^| than the basis of 
the internal problem used up to now (see also [36^. The 
latter matrix is diagonalized by an orthogonal transfor- 



mation rj 



1 77W 77(2) 



.77^), where each entry 17 w is a 
real vector in the L-dimensional part of the total Hilbert 
space, which corresponds to the considered bunch. Due 
to the factorized structure, the interaction matrix pos- 
sesses the only nonzero eigenvalue T c m ^ = A 2 , which 
belongs to the eigenvector tjW = a = A/V A 2 . As a 
result, the effective Hamiltonian obtains in the doorway 
basis the form (we drop below the index m of the bunch) 



n eff — 



(102) 



h e 

The following notation has been used here (/x=2,3,...,L): 



^E^ (D W-2A L (m^) 2 H4 D) > 



5e r 



(103) 



with the L—l vectors rj^) and elements e u of the diagonal 
matrix e being defined by the eigenvalue problem 



E f (») „W = e S 
fc n 'In 'In ~ t fJ- u V-v 



(104) 



in the (L— l)-dimensional subspace orthogonal to the vec- 
tor a. The quantity (e[ D ^ ) is the weighted mean position 
of the internal levels. The non-Hcrmitian Hamiltonian 
(|102f> describes a wide doorway resonance with width T c 
displaced from the bunch by the distance 8e c and coupled 
to a background of L — 1 stable states by means of the 
matrix elements h n . Only due to this interaction do such 
states get access to the continuum via the doorway state 
existing in the outer well. 

With the help of the completeness condition 



EM (/-») 



at ai> 



one finds from Eq. (|1U4|I 



'h ~ "7< 



(11 . 



The orthogonality condition a ■ r\ 
gives the equation 



-(D) 



(105) 



(106) 



immediately 



(107) 



for the new positions of the stable levels. This equation 
shows that each new level e a lies between two neighboring 

old ones e| D ^ and therefore is shifted with respect to the 
latter only by a distance comparable with the initial mean 
level spacing. This is much smaller than the displacement 
Se c of the collective level e c . 

The interaction h mixes these states and forms L fi- 
nal resonance states. The complex energies of the exact 
states, i.e. the eigenvalues of the effective Hamiltonian 
l)102[l . satisfy the secular equation 



£ 



1 

- r 

2 u 1 



h 2 



(108) 



In particular, for the collective doorway state one obtains 
from this equation 



£coll 



-(D) v 



,2 1 



1 - 



(Se c ) 2 



Se c - 



,2 1 



1 - 



(Se c ) 2 



r c 

(109) 



Equations (|103(l and i|105fl allow us to express the square 
length h 2 of the mixing vector in terms of the variance 
of the non-perturbed internal levels, 

h 2 = B^af-fcfV) 



= [Ae (D) ] 2 . 



(110) 



Thus, interaction is weak in the doorway basis under the 
condition 8e c ^> Ae' D ' and the doorway state keeps al- 
most the whole energy displacement and width. All other 
states are trapped in the interior region and share small 
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portions ~ [Ae' D '/i5e c ] 2 of the collective displacement 
and width in accordance with the group velocity |23j . 
The solution just described has formally very much in 
common with the schematic model of the so-called nu- 
clear dipole giant resonance developed in [3(| |!3 • 



V. CONCLUSION 

In this paper we analyzed the relevance of the con- 
cept of the non-Hermitian effective Hamiltonian in finite- 
range potential scattering. Single-channel s-wave reso- 
nance scattering is considered as an example. The scat- 
tering function S(k) is meromorphic in this case, i.e., has 
only isolated poles in the complex plane of the wave num- 
bers. The number of poles is, generally, infinite but only 
a finite part of them can be interpreted as resonances. 

We first presented (Sees. II and III) a consistent self- 
adjoint formulation of the scattering problem, which is 
based on separation of the configuration space into inter- 
nal and external segments. In this way, the S(k) func- 
tion is represented in terms of a non-Hermitian energy- 
dependent operator H(k). There exists a wide freedom in 
choosing the radius a of separation as well as the bound- 
ary conditions at this point. Different choices yield dif- 
ferent explicit forms of this operator together with the 
S(k) and R(k) functions. This can, in particular, come 
out strongly if one truncates the matrix Ti to calculate 
numerically the positions of the poles of S function in the 
complex k plane. Nevertheless, the true complex poles of 
the S function and this function itself depend, as we ex- 
plicitly demonstrate, neither on the BC nor on the radius 
a. 

Although all choices are formally allowed, this does 
not mean that they are equally adequate from the phys- 
ical point of view. For instance, a fictitious immediate 
reflection takes place at the (arbitrarily chosen) point 
of separation a. The artificial separation of the phase 
of this reflection results, in turn, in the appearance of 
an infinite number of remote poles of the function R(k), 
which describe the smooth contribution from the inte- 
rior region r < a. The phase of the immediate reflection 
must be compensated by the contributions of such poles 
for the radius a to disappear finally from the scattering 
amplitude. This can cause unjustified complications in 
intermediate stages of calculations. We argue that the 
harm is minimized if the separation radius matches an 
outer potential barrier when it is strong enough to make 
immediate reflection at this point quite probable. 

The notion of the non-Hermitian effective Hamiltonian 
first introduced in the theory of resonance nuclear re- 
actions emerges when a group of very close resonance 
states strongly overlap and interfere. We stress that, in 
contrast to the nuclear reactions, in the cases of the po- 
tential resonance scattering usually discussed in the lit- 
erature the density of the energy spectrum is too low 
and the ordinary Breit-Wigner approximation of isolated 
resonances usually suffices. The complex energy of the 



roth resonance which dominates in the scattering ampli- 
tudes near the scattering energy E w \z m \ 2 is well ap- 
proximated by the diagonal matrix element T~L mm {\z m \). 
There is no room for an effective Hamiltonian in this 
case or, more strictly, it is embodied by the lxl "ma- 
trix" Timrn- We emphasize that the energy dependent 
operator H(k) should not be confused with the effective 
Hamiltonian. Energy dependent eigenvalues of this op- 
erator are not, generally, in one-to-one correspondence 
with the complex energies of the actual resonance states. 

Overlap and interference of the resonances can become 
possible in the cases when the energy spectrum has a 
band structure. As an example of this kind we investi- 
gate a periodically disposed chain of a finite number L 
of radial S barriers. All 5 1 , K and R functions are found 
in closed forms in terms of an L x L fc-dependent matrix 
propagator. This allows us to study in detail all ana- 
lytical properties in the complex k plane and to verify 
the correspondence with the projection formalism used. 
There exist a finite number of separated bands of close 
resonances. Within the roth band which lie near the 
scattering energy E w \z m \ 2 one can neglect all smooth 
variations with k. In this approximation, the mentioned 
propagator proved to coincide with the resolvent of the 
L x L block TLii'(\z m \) of the matrix TimL(k) taken at 
the fixed value k = \z m \. Just this matrix plays the role 
of the effective Hamiltonian of the system in the energy 
interval within the band. 

Different choices of BCs yield different patterns of the 
resonance interference. In particular, the spectrum of 
the poles of the R function is exactly reproduced in the 
framework of the projection technique with the BC l|51|) 
of Neumann type fixed at the position a = L of the outer 
barrier. The corresponding energy levels depend on the 
penetrability constant A L of the outer barrier. Shifts of 
the levels due to the coupling to the continuum are in- 
cluded in this case from the very beginning. 

Utilization of the Dirichlet BC for the intrinsic motion 
gives another but yet equivalent formulation of the scat- 
tering problem considered. The internal problem in this 
case fixes a closed counterpart of the open system under 
consideration. This enables us to investigate the change 
of the regime of the internal motion of the system as its 
openness grows. The interaction via the continuum shifts 
the original A L -independent levels of the internal motion 
along both real and imaginary axes. A transition is ex- 
plicitly demonstrated from the bunch of L similar narrow 
resonances to the formation of a relatively broad reso- 
nance strongly shifted with respect to the band, which 
exists in the outer well. The other resonances turn out to 
be trapped in the inner part of the system. This is quite 
similar to the results of earlier investigations which rely 
upon the notion of the non-Hermitian effective Hamilto- 
nian of an open system. 



16 



Acknowledgments 

We are grateful to Y.V. Fyodorov for discussions on 
Feshbach's projector technique in an early stage of this 
work, to H. Schanz and V. Zelevinsky for interesting dis- 
cussions, and to G. Hackenbroich and C. Viviescas for 
their interest in the work and useful conversations. One 



of us (V.V.S.) greatly appreciates the generous hospital- 
ity of the Max-Planck Institute for the Complex Systems 
extended to him during his stay in Dresden. The financial 
support by RFBR Grant No. 03-02-16151 (D.V.S. and 
V.V.S.) and SFB 237 "Unordnung and grosse Fluktuatio- 
nen" (D.V.S. and H.J.S.) is acknowledged with thanks. 



[1] H. Feshbach, Ann. Phys. (N.Y.) 5, 357 (1958). 

[2] H. Feshbach, Ann. Phys. (N.Y.) 19, 287 (1962), 
[reprinted in Ann. Phys. (N.Y.) 281, 519 (2000)]. 

[3] P. L. Kapur and R. Peierls, Proc. R. Soc. London, Ser. 
A 166, 277 (1938). 

[4] E. Wigner and L. Eisenbud, Phys. Rev. 72, 29 (1947). 

[5] M. S. Livsic, Sov. Phys. JETP 4, 91 (1957). 

[6] M. S. Livsic, Operators, Oscillations, Waves (Open Sys- 
tems), Translations of Mathematical Monographs Vol. 34 
(American Mathematical Society, Providence, RI, 1973). 

[7] C. Mahaux and H. A. Weidenmiiller, Shell-Model Ap- 
proach to Nuclear Reactions (North-Holland, Amster- 
dam, 1969). 

[8] J. J. M. Verbaarschot, H. A. Weidenmiiller, and M. R. 

Zirnbauer, Phys. Rep. 129, 367 (1985). 
[9] V. V. Sokolov and V. G. Zelevinsky, Nucl. Phys. A 504, 

562 (1989). 

[10] F. Remacle, M. Munster, V. B. Pavlov- Verevkin, and 
M. Desouter-Lecomte, Phys. Lett. A 145, 265 (1990); M. 
Desouter-Lecomte and F. Remacle, Chem. Phys. 164, 11 
(1992). 

[11] I. Rotter, Rep. Prog. Phys. 54, 635 (1991). 

[12] Y. V. Fyodorov and H.-J. Sommers, J. Math. Phys. 38, 

1918 (1997). 
[13] F.-M. Dittes, Phys. Rep. 339, 215 (2000). 
[14] Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 
[15] H.-J. Stockmann, Quantum Chaos: An Introduction 

(Cambridge University Press, Cambridge, UK, 1999). 
[16] G. Akguc and L. E. Reichl, Phys. Rev. E 64, 056221 

(2001). 

[17] K. Pichugin, H. Schanz, and P. Seba, Phys. Rev. E 64, 
056227 (2001). 

[18] H.-J. Stockmann, E. Persson, Y.-H. Kim, M. Barth, 
U. Kuhl, and I. Rotter, Phys. Rev. E 65, 066211 (2002). 

[19] G. Hackenbroich, C. Viviescas, and F. Haake, Phys. Rev. 
Lett. 89, 083902 (2002); C. Viviescas and G. Hacken- 
broich, Phys. Rev. A 67, 013805 (2003). 

[20] H. Alt, P. von Brentano, H.-D. Graf, R.-D. Herzberg, 
M. Philipp, A. Richter, and P. Schardt, Nucl. Phys. A 
560, 293 (1993). 

[21] H. Alt, H.-D. Graf, H. L. Harney, R. Hofferbert, 
H. Lengeler, A. Richter, P. Schardt, and H. A. Wei- 
denmiiller, Phys. Rev. Lett. 74, 62 (1995). 

[22] E. Persson, I. Rotter, H.-J. Stockmann, and M. Barth, 



Phys. Rev. Lett. 85, 2478 (2000). 
[23] V. V. Sokolov and V. G. Zelevinsky, Ann. Phys. (N.Y.) 

216, 323 (1992). 
[24] T. Kottos and U. Smilansky, Ann. Phys. (N.Y) 274, 76 

(1999). 

[25] F. Barra and P. Gaspard, J. Phys. A: Math. Gen. 32, 

3357 (1999). We are grateful to T. Kottos for drawing 

our attention to this paper. 
[26] C. Bloch, Nucl. Phys. 4, 503 (1957). 
[27] A. M. Lane and D. Robson, Phys. Rev. 151, 774 (1966); 

D. Robson and A. M. Lane, ibid. 161, 982 (1967). 
[28] A. M. Lane and R. G. Thomas, Rev. Mod. Phys. 30, 257 

(1958). 

[29] As usual, factors e~ an with infinitesimal a should be 
introduced to guarantee the absolute convergence. 

[30] Let us recall that the penetrability of a S— barrier D(k) = 
(Afc) 2 /[1 + (Afc) 2 ] is small, D(k) » (Afc) 2 , if Afc < 1, and 
goes to one under the opposite condition. 

[31] N. Hu, Phys. Rev. 74, 131 (1948). 

[32] H. M. Nussenzveig, Causality and Dispersion Relations 

(Academic Press, New York, 1972). 
[33] E. Persson, M. Miiller, and I. Rotter, Phys. Rev. C 53, 

3002 (1996). 

[34] V. Zelevinsky and A. Volya, in Challenges of Nuclear 
Structure, edited by A. Covello (World Scientific, Singa- 
pore, 2002), p. 261; see also A. Volya and V. Zelevinsky, 
e- print nucl-th/02 11039 

[35] By a slight shift of the position of the outer barrier 
and subsequent careful tuning of the two free parameters 
complete degeneracy of the resonances can be achieved; 
see J. S. Bell and C. J. Goebel, Phys. Rev. 138, B 1198 
(1965); E. Hernandez, A. Jauregui, and A. Mondragon, 
J. Phys. A: Math. Gen. 33, 4507 (2000). Recently, the 
experimental study of such an exceptional point has been 
performed in a microwave cavity setup: C. Dembowski, 
H.-D. Graf, H. L. Harney, A. Heine, W. D. Heiss, H. Re- 
hfeld, and A. Richter, Phys. Rev. Lett. 86, 787 (2001); 
see, however, I. Rotter, Phys. Rev. E 67, 026204 (2003). 

[36] V. V. Sokolov, I. Rotter, D. V. Savin, and M. Miiller, 
Phys. Rev. C 56, 1031 (1997). 

[37] V. V. Sokolov and V. G. Zelevinsky, Fizika (Zagreb) 22, 
303 (1990). 



